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ABSTRACT 

We present an investigation of satellite galaxies in the outskirts of galaxy clusters 
taken from a series of high- resolution iV-body simulations. We focus on the so-called 
"backsplash population" , i.e. satellite galaxies that once were inside the virial radius 
of the host but now reside beyond it. We find that this population is significant in 
number and needs to be appreciated when interpreting the various galaxy morphol- 
ogy environmental relationships and decoupling the degeneracy between nature and 
nurture. Specifically, we find that approximately half of the galaxies with current clus- 
tercentric distance in the interval 1-2 virial radii of the host are backsplash galaxies 
which once penetrated deep into the cluster potential, with 90% of these entering to 
within 50% of the virial radius. These galaxies have undergone significant tidal dis- 
ruption, loosing on average 40% of their mass. This results in a mass function for the 
backsplash population different to those galaxies infalling for the first time. We further 
show that these two populations are kinematically distinct and should be observable 
within existent spectroscopic surveys. 
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1 INTRODUCTION 

The relationship between galaxy morphology and local en- 
vironment (i.e. the morphology-density relation) was first 
noticed by Hubble & Humason (1931), where they re- 
ported that field and cluster galaxy populations differ. Oem- 
ler (1974) extended this finding by showing that the rela- 
tionship held for differing clusters richness. The field truly 
emerged when Dressier (1980) demonstrated the strong re- 
lationship over five orders of magnitude between the local 
density of galaxies and the proportions of different morpho- 
logical types. Bhavsar (1981), de Souza et al. (1982), and 
Postman & Geller (1984) extended this work further to in- 
clude the environments of both loose and compact groups. 
Recently, Aguerri et al. (2004) performed a thorough anal- 
ysis of 116 bright galaxies in the Coma cluster, finding that 
bluer objects are located at larger projected radii while si- 
multaneously showing a larger velocity dispersion than their 
red counterparts. Moreover, the bluest objects also host the 
most prominent disks contrary to systems observed close 
to the cluster centre or in high-density environments. Envi- 
ronmental dependence of galactic stellar populations is also 
seen in the Butcher-Oemler effect (Butcher & Oemler 1978, 
Kodama & Bower 2001) with clusters at higher redshifts 
showing a greater fraction of blue objects than are seen at 
present. 



The above observational work supports the idea that 
galaxies in clusters are substantially different from galaxies 
in the field. But the origins of these morphology-density re- 
lationships are still not fully understood with several large 
and small scale mechanisms proposed to explain their exis- 
tence, including ram pressure stripping (Gun & Gott 1972), 
tidal stripping/star formation (Merritt 1983, 1984), starva- 
tion (Larson, Tinsley & Caldwell 1980), galaxy merger and 
harassment (Icke 1985; Moore et al. 1996) and dynamics 
(Tsuchiya & Shimada 2000). 

Recent cosmological simulations (both hydrodynami- 
cal and iV-body) have provided valuable insight into the 
mechanisms responsible for the morphology-density rela- 
tionship (Springel et al. 2001; Goto et al. 2003; Okamoto 
& Nagashima 2003). In the analysis which follows, we fo- 
cus on the dynamics of satellite galaxies taken from a 
series of high-resolution, fully self-consistent, cosmological 
simulations of eight galaxy clusters. We concentrate on 
the outskirts of these clusters, i.e. distances in the range 
[R v ir,2.5R V ir], which (observationally) have only recently 
been probed through wide-field optical imaging and spec- 
troscopy (Miyazaki et al. 2002; Lewis et al. 2002). We will 
demonstrate that a rich population of galaxies exist beyond 
the virial radius most of which have previously spent time 
near the cluster centre and can be seen in Figure 1. We 
characterise the spatial, velocity and mass properties of this 
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Figure 1. A simulated cluster at z=0 with the virial radius indi- 
cated by the dark sphere. The other line in this figure represents 
the orbital path of a "backsplash" galaxy. Such a galaxy has pre- 
viously spent time near the cluster centre but now lies outside 
the virial radius of the cluster. 

population and contrast these with those of the spatially 
coincident newly infalling galaxies. 

Our work complements the earlier studies of Balogh 
et al. (2000) and Mamon ct al. (2004). Balogh et al. investi- 
gated the particle backsplash from cosmological simulations 
and found that 50 ± 20% of the particles within [R200, 2.R200] 
had passed through the -R200 radius. Mamon et al. recently 
extended this work to calculate the maximum backsplash 
distance for particles to be 2.5-Rioo and for galaxies 1.7-Rioo- 

The outline of the paper is as follows. Section 2 provides 
a description of the cosmological simulations employed. We 
investigate the number distribution of galaxies in the cluster 
outskirts in Section 3, with the mass distribution discussed 
in Section 4 and the velocity distribution of the satellites 
investigated in Section 5. We finish with our summary and 
conclusions in Section 6. 



2 SIMULATION DETAILS 

Our analysis is based on a suite of eight high-resolution N- 
body simulations (Gill, Knebe & Gibson 2004a) carried out 
using the publicly available adaptive mesh refinement code 
MLAPM (Knebe, Green & Binney 2001) in a standard ACDM 
cosmology (Q = 0.3, Q x = 0.7, Q. b h 2 = 0.04, h = 0.7, a 8 = 
0.9). Each run focuses on the formation and evolution of a 
dark matter galaxy cluster containing of order one million 
particles, with mass resolution 1.6 x 10 s ft -1 Mq and force 
resolution ~2/i _1 kpc which is of the order 0.05% of the 
host's virial radius. These simulations have the required res- 
olution to follow the satellites within the very central regions 
of the host potential (>5-10% of the virial radius) and the 
time resolution to resolve the satellite dynamics with good 



accuracy (At «170 Myrs). Such temporal resolution pro- 
vides of order 10-20 timesteps per orbit per satellite galaxy, 
Thus allowing these simulations to be used in a previous 
paper Gill et al. (2004b) to accurately measure the orbital 
parameters of each individual satellite galaxy. 

The clusters were chosen to sample a variety of environ- 
ments and their details are summarised in (Gill, Knebe & 
Gibson 2004a). We define the virial radius i? v i r as the point 
where the density of the host (measured in terms of the 
cosmological background density pi,) drops below the virial 
overdensity A v i r = 340. A v ; r = 340 is based upon the dissi- 
pationless spherical top-hat collapse model and is a function 
of both cosmological model and time. We further applied a 
lower mass cut for all the satellite galaxies used in this pa- 
per at 2 x 10 10 h^ 1 Mq (100 particles). For further specific 
details of the host halos such as masses and density pro- 
files please refer to the earlier papers in the series (cf. Gill, 
Knebe & Gibson 2004a and Gill et al. 2004b). 

For later reference and comparison to previous studies 
by Balogh et al. (2000) and Mamon et al. (2004) we relate 
R v ir to other definitions of the virial radius, namely -R200 
and i?ioo- According to the NFW profile (Navarro, Frenck & 
White 1997) the dark matter density in the outer regions of 
halos drops like p oc r -3 and hence R vir ~ I.4.R200 and 

Rvir ~ -Rl00- 



3 NUMBER DISTRIBUTION OF GALAXIES 
IN THE CLUSTER OUTSKIRTS 

The fact that galaxies in clusters are different from galaxies 
in the field has been attributed to a variety of mechanisms 
acting on a variety of scales. It has been suggested that 
galaxies can "rebound" up to 1.7J?ioo (Mamon et al. 2004), 
thus blurring the definition of "cluster" and "field". In 
this section we investigate cluster galaxy "backsplash" : how 
significant is the population of galaxies that once passed 
through a cluster's virial radius but now resides in its out- 
skirts? 

In Figure 2 we plot for all of our eight simulated clusters 
the minimum distance D m in a galaxy reached to the clus- 
ter centre throughout its history versus its current distance 
D z= o- Both distances have been normalised by the cluster's 
present day virial radius i? v ir- There are four distinct popu- 
lations of satellites visible in this figure: 

1. Pmin = Dz=o'- (the infalling population) 
These satellites are falling in for the first time. 

2. Dmin > Rvir and D z= p > -Rvj r : (infalling sub-population) 
This population is made up of subhalos orbiting satellites, 
i.e. sub-subhalos. 

3. Dmin < Rvii and Dz=o > Rvir'- (backsplash satellites) 
These satellites once passed through the virial radius of the 
respective host but have "rebound" to the outskirts of the 
galaxy cluster. 

4. Dmin < Rvii and Dz=o < Rvir'- (the bound population) 
This is the "normal" satellite population orbiting within the 
virial radius and bound to the host. 
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Figure 2. Minimum distance for all satellite galaxies versus current distance in terms of the host's virial radius. The population in the 
upper right corner consists of substructure orbiting within satellites, i.e. sub-subhalos. This figure clearly indicates that there is a distinct 
population of "backsplash" satellites, i.e. D m ; n < R y ^ I and .D z= o > i? v ir m the lower right corner. 
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Figure 3. The relative normalised distribution function of the 
minimum distance D mln /R v i T for the whole galaxy popula- 
tion (thick histograms) and the satellites with current positions 
greater than the virial radius (thin histograms). We note that 
they both peak near 25% of the cluster virial radius. 



Our data sets indicates that the maximum backsplash 
distance encountered is ~ 2.5R V i r which is larger than the 
value presented in Mamon et al. (2004)* yet not widely in 
disagreement, further very few backsplash galaxies exist be- 
yond ~ 2R v i r . 

A more quantitative analysis reveals that 30% of all 
satellites with -D m in < -Rvir are outside the host's virial ra- 
dius at z = 0. Further, 50% of all galaxies with a current 
distance D z= o in the range [.Rvir, 2-R v i r ] are in fact back- 
splash galaxies, consistent with the value 50 ± 20 % quoted 
by Balogh et al. (2000) (based on cluster particles as opposed 
to gravitationally bound satellites though). 

Only 2% of the backsplash population has had more 
than one orbit (i.e. eight satellites in total) and most stem 
from single passages through the host's virial radius, respec- 
tively. And each of these eight "multiple passage" backsplash 
galaxies are found in the range [-Rvir, 1.2R v i r ] which is close 
to the host. 

We now investigate the depth to which backsplash 
galaxies penetrate the cluster potential. To this extent we 
plot in Figure 3 the normalised number distribution of the 
minimum distance D m in for the entire satellite population 
(thick histogram) and the restricted set of satellites cur- 
rently outside the virial radius (thin histograms). It is in- 

* refer to Section 2 for conversion of R v ir to Rioo- 
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Figure 4. Relative mass loss for the backsplash satellites over the 
lifetime of the host as a function of minimum distance. The solid 
line shows the average with the dashed lines being the standard 
deviation. 



Figure 5. The relative number distribution of satellite masses. 
The thin histogram represents the backsplash population and the 
thick histogram shows the satellites infalling for the first time. 



teresting to note that the distributions are fairly similar for 
£>mm < Rvir, i.e. they both peak near 0.25. We stress that 
this is the minimum distance as measured over the lifetime 
of the host halo and hence is not to be confused with the 
last pericentre whose distribution peaks at about 35% of 
i?vir (Gill et al. 2004b). Figure 3 reveals that the a number 
of the galaxies that are now in the outskirts of the clus- 
ter have had passages as close to the host as their bound 
counterparts which implies that the former are on highly 
radial orbits. This is consistent with the findings of Solanes 
et al. (2001) who found evidence that gas-poor spirals in HI 
deficient clusters move in orbits more radial than those of 
the gas-rich objects. In fact, 90% of our backsplash satellites 
once passed within the inner 50% of the virial radius. 

In summary the number of backsplash galaxies is sig- 
nificant and should be accounted for when interpreting the 
galaxy morphology-density relationship. These galaxies pen- 
etrate deep within the cluster potential, as deep as their 
bound counterparts. Hence, they should be sampling the 
large and small-scale transformation mechanisms alluded to 
in Section 1. For example, if we used the prescriptions out- 
lined in Treu et al. (2003), these galaxies would also have un- 
dergone starvation, ram pressure stripping, tidally triggered 
star formation and significant tidal stripping, thus confus- 
ing any correlation of the galaxy's morphology with local 
environment. 



4 MASS DISTRIBUTION OF THE GALAXIES 
IN THE CLUSTER OUTSKIRTS 

One way to gauge the significance of these transformation 
mechanisms is to investigate the mass loss of a galaxy during 
its encounter with the cluster. To this extent we focus only 
on the backsplash population and plot the relative mass lost 
over the lifetime of the host as a function of minimum dis- 
tance -D m in/-Rvir for each of these satellite in Figure 4. We 



find that the closer a satellite gets to the host the stronger 
the tidal stripping. This result is not surprising but needs 
to be viewed from the perspective of the backsplash popu- 
lation: galaxies in the outskirts experienced significant tidal 
interactions and mass loss, in some cases up to 80% of their 
original mass. 

Further to that we see little dependence of present 
day clustercentric distance on the mass lost (not presented 
though). At all distances outside the virial radius the aver- 
age mass lost for each backsplash galaxy is ~ 40%. 

Since the population of backsplash galaxies has under- 
gone significant tidal stripping we expect the mass spectrum 
to be different from that of the galaxies infalling for the first 
time. In Figure 5 we plot the relative number distribution of 
satellites with a certain mass M aat measured in terms of the 
host's virial mass M host . The thin histogram represents the 
backsplash population and the thick histogram shows the in- 
falling population. To better discriminate between these two 
populations we fitted the mass spectra to a simple power-law 

n{M) = CM a . (1) 

There is a marginal difference between the slopes, the back- 
splash population having a — 0.9 ±0.3 whereas the infalling 
satellites distribution is characterised by a = 0.7 ± 0.3. The 
steeper slope of the backsplashed galaxies reflects the lack 
of larger satellites and an over-abundance of smaller galax- 
ies, respectively. This tilt in the mass spectrum is readily 
explained by the mass loss highlighted in Figure 4. We note 
that the first mass bin was not used in fitting the power- 
law as it simply reflects the absolute mass cut of 2 x 10 10 
h~ x M (i.e. we did not use a relative mass cut with respect 
to the mass of the host). 
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5 VELOCITY PROPERTIES OF GALAXIES IN 
THE CLUSTER OUTSKIRTS 

To populate the outskirts of clusters, galaxies that once 
passed through the virial radius and close to the centre of 
the host must have had high velocities. In this section we 
explore the relative velocities of satellite galaxies within the 
cluster and its outskirts. 



Table 1. Number of inbound and outbound backsplash galaxies 
in different radial bins. 



distance receding galaxies approaching galaxies 



(1.00 - 1.33) R vir 102 106 

(1.33 - 1.67) R viT 87 57 

(1.67 - 2.00) R vil 47 22 



5.1 A kinematically distinct backsplash 
population 

Figure 6 shows the absolute value of the relative velocity 
between a satellite galaxy and the host cluster as a function 
of current clustercentric distance. We present the results for 
two populations again, namely the infalling satellites that 
are still outside the virial radius (diamonds) and the back- 
splash galaxies alongside galaxies within the virial radius 
(crosses) . 

At the virial radius the maximum relative velocity for a 
galaxy in our data set is 1500 km s _1 with no galaxy farther 
away having a greater velocity. This finding can be used 
to explore the two competing models for the ram-pressure 
stripping observed in NGC 4522, a spiral galaxy in the Virgo 
cluster with a distance close to the virial radius (Vollmer 
ct al. 2001). Kenney et al. (2004) suggested two possible 
scenarios to explain its ongoing ICM-ISM stripping: either 
NGC 4522 is experiencing stripping at this distance due to 
bulk motions and local density enhancements of the ICM 
produced by shocks in the ICM, or - since ram-pressure 
stripping is proportional to v 2 - the galaxy has a very high 
relative velocity, of order 4000 km s _1 . As hinted by Kenney 
ct al., such a velocity is extremely unlikely, and as shown 
in Figure 6 only in the very central regions of clusters do 
galaxies reach these high relative velocities. 

Figure 6 furthermore shows that the backsplash popu- 
lation is quite distinct kinematically; the infalling satellites 
have significantly larger velocities than the backsplash galax- 
ies. This provides a possible mechanism to observationally 
detect these rebound satellites. 

To investigate this in more detail we divided all (in- 
falling and backsplash) satellites into three radial bins be- 
tween 1 and 2 R v i T . In Figure 7 we show the distribution of 
the relative velocities for all satellites in the respective bin 
(solid lines) . The dotted (dashed) lines are for the backsplash 
(infalling) population alone, normalised to the total number 
of satellites (i.e. the actual distribution is simply the sum of 
these two distributions). From the distributions presented 
in Figure 7 it is rather obvious that the backsplash popu- 
lation should be detectable by simply plotting the velocity 
distribution function (VDF) for cluster galaxies in the range 
[i? v ir, 2ii v i r ]: there are two distinct peaks with the lower ve- 
locity peak indicative of the rebound satellites. Moreover, as 
we move further away from the host the separation between 
the two populations becomes even more pronounced. 

One might still pose the question though, if the back- 
splash population can be identified by their "flight path" , i.e. 
the orientation of their velocity vector. In Table 5.1 we there- 
fore summarise the numbers of backsplash galaxies moving 
towards ("approaching") and away ("receding") from the 
centre of the host. This demonstrates that the backsplash 
galaxies cannot simply be selected by the direction of their 
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Figure 7. The distribution of the relative velocities in three dis- 
tance bins, i.e. 1.00-1.33, 1.33-1.67 and 1.67-2.00 virial radii. The 
dotted histograms shows the backsplash galaxies while the dashed 
lines represent the infalling galaxies (both normalized to the to- 
tal number of satellites). The solid histogram is the sum of both 
distributions. We see that the distribution has two distinct peaks. 



velocity as they appear to move in both directions. How- 
ever, as distance increases we do find a higher percentage of 
backsplash galaxies moving away from the host. 



5.2 Observational impact 

Thus far, we have relied upon the fact that our simulation 
data provides full six dimensional velocity and spatial infor- 
mation. In this section we extend the velocity distribution 
function to the "observer's plane" by restricting the data to 
projected distances and line-of-sight velocities; we make this 
transformation by placing the potential observer at infinity. 

In Figure 8 we show the distribution of the line-of-sight 
velocities (still with respect to the host) for galaxies with 
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Figure 6. The absolute value of the relative velocity for the infalling satellites (diamonds) and the population of satellite that have 
crossed the virial radius at sometime (crosses) as a function of current cluster distance to the host. The middle solid lines represent 
averages with the outer lines being the standard deviations. Note the kinematic distinction between the populations. 



projected distances between 1.0 and 2.0 virial radii. Veloci- 
ties have also been convolved with a velocity uncertainty of 
100 km s _1 , typical of that encountered in existing multi- 
object spectroscopic surveys such as the 2dFGRS (Colless 
et al. 2001). The dotted histogram shows the backsplash 
galaxies while the dashed histogram represents the infalling 
galaxies (both normalised to the total number of satellites). 
The solid histogram is the sum of both distributions. 

Not surprisingly, this figure appears quite different from 
Figure 7 due to the loss of information when transforming 
from 6-dimensional space to 3-dimensional space. The kine- 
matically distinct bimodal populations are no longer read- 
ily separable. However, if the backsplash population does 
not exist, one should simply observe the dashed histogram 
in Figure 8. The presence of a backsplash population "dis- 
torts" the Gaussian velocity distribution function such that 
it increases with decreasing line-of-sight velocity. 

Figure 8 should now be comparable to extant cluster 
data sets such as the 2dFGRS (Lewis et al. 2002). Lewis 
et al. show VDFs for 17 clusters (of order 50 galaxies per 
cluster). Stacking the 17 cluster VDFs, eliminating all galax- 
ies within the (projected) virial radius, should yield a com- 
posite 1-2 R vir VDF with ~ 200 — 300 galaxies; such a 
composite VDF should either support or refute our pre- 
dicted result (Figure 8). However, this is just a first step 



to a fair comparison. The next step would require colour, 
surface brightness and kinematic selection. 



6 CONCLUSIONS 

Observational data supports the idea that galaxies in clus- 
ters are substantially different from galaxies in the field. 
There is a clear correlation between galaxy morphology and 
density of the local environment. However, the origin of this 
relation is far from being understood. 

In this study we have presented an analysis of satellite 
galaxies that once passed through the virial radius close to 
the centre of their respective host halo, but are now found 
outside the virial radius in the outskirts of the cluster. We 
have shown that this backsplash population is not negligible 
and needs to be accounted for when interpreting the various 
galaxy morphology relationships and decoupling the degen- 
eracy between nature and nurture. 

We must also appreciate that the infalling population 
is not expected to be pristine. Rather, we would expect that 
infall galaxies have undergone some sort of pre-processing 
in groups before entering the cluster too as indicated by the 
sub-subhalos in Figure 2. 

Our results can be summarized as follows: 

• 30% of all galaxies that ever came closer to the 
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Figure 8. The distribution of the line-of-sight velocities (relative 
to the host) for galaxies with projected distance between 1.0-2.0 
virial radii. Velocities have been convolved with the 2dF velocity 
uncertainty of 100 km s _1 , typical of that encountered in multi- 
object spectroscopic surveys (e.g. Colless et al. 2001). The dotted 
histogram represents the backsplash galaxies while the dashed 
histogram is based upon the infalling galaxies (both normalized 
to the total number of satellites). The solid histogram is the sum 
of both distributions. 



host than its virial radius are now located in the range 

vir ■ 2.5R vir I , 

• 50% of all galaxies in the region [J? v i r ,27i v j r ] are back- 
splash galaxies, 

• 90% of the backsplash galaxies penetrated deeper than 
50% of R V [ T into the host's potential, 

• during their passage through the cluster, on average the 
backsplash galaxies lose 40% of their mass, thus 

• the mass spectrum of the backsplash population has 
a steeper power-law slope than their infalling counterparts, 
thus it has fewer massive galaxies and more light ones, 

• the velocities of the infalling satellites is too small to 
account for ram-pressure stripping in the cluster outskirts, 

• the backsplash population has a factor of two smaller 
relative velocity than the infalling satellites, making it kine- 
matically distinct. 

When transforming the last result into the observers 
plane though, the velocity separation between the infalling 
and backsplash population is removed. However, the back- 
splash population should still be detectable as it is responsi- 
ble for a continuous rise in the distribution function towards 
low line-of-sight velocities. 

Our results suggest that we not only expect the back- 
splash population to experience various large-scale transfor- 
mation mechanisms, but also small-scale ones, undergoing 
starvation, ram pressure stripping, tidally triggered star for- 
mation and significant tidal stripping. 
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